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Structures of calcium peroxide (CaCh) are investigated in the pressure range 0-200 GPa using 
the ab initio random structure searching (AIRSS) method and density functional theory (DFT) 
calculations. At 0 GPa, there are several CaC >2 structures very close in enthalpy, with the ground- 
state structure dependent on the choice of exchange-correlation functional. Further stable structures 
for CaC >2 with C2/c, /4/mcm and P2\/c symmetries emerge at pressures below 40 GPa. These 
phases are thermodynamically stable against decomposition into CaO and O 2 . The stability of Ca02 
with respect to decomposition increases with pressure, with peak stability occurring at the CaO Bl- 
B2 phase transition at 65 GPa. Phonon calculations using the quasiharmonic approximation show 
that Ca 02 is a stable oxide of calcium at mantle temperatures and pressures, highlighting a possible 
role for Ca02 in planetary geochemistry. We sketch the phase diagram for Ca02, and find at least 
five new stable phases in the pressure/temperature ranges 0 < P < 60 GPa, 0 < T < 600 K, 
including two new candidates for the zero-pressure ground state structure. 

PACS numbers: 


INTRODUCTION 

The typical oxide formed by calcium metal is calcium 
oxide, CaO, having Ca and O in +2 and -2 oxidation 
states respectively. Calcium and oxygen can also com¬ 
bine to form calcium peroxide, Ca 02 , a compound which 
enjoys a variety of uses in industry and agriculture. Cal¬ 
cium peroxide is used as a source of chemically bound 
but easily evolved oxygen in fertilisers, for oxygenation 
and disinfection of water, and in soil remediation [1, 2]. 

At ambient pressure bulk calcium peroxide decom¬ 
poses at a temperature of about 620 K [2, 3]. Early 
X-ray diffraction (XRD) experiments assigned a tetrago¬ 
nal ‘calcium carbide’ structure of space group 14,/mmm 
to CaC >2 [4]. This same structure was already known 
to be formed by heavier alkaline earth metal peroxides 
[5, 6]. Recently, Zhao et al. [7] used an adaptive genetic 
algorithm and density functional theory (DFT) calcula¬ 
tions to search for structures of Ca 02 , finding a new or¬ 
thorhombic ground state structure of Pna2\ symmetry, 
which is calculated to be close to thermodynamic sta¬ 
bility at zero pressure and temperature. The simulated 
XRD pattern from this structure is in good agreement 
with the available experimental data [7]. Thermody¬ 
namic stability in this case means stability against de¬ 
composition via 

Ca0 2 —> CaO + io 2 . (1) 

We are interested in the stabilities of structures of 
Ca02 from zero pressure up to 200 GPa, and temper¬ 
atures up to 1000 K. Calcium and oxygen have high 
abundances in the Earth’s crust and mantle and, be¬ 
cause they also have high cosmic abundances, stable com¬ 


pounds formed from these elements at high pressures are 
key (exo)-planetary building blocks. Understanding the 
structures of such compounds allows insight into the com¬ 
position of planetary interiors, including exoplanets. To 
date, almost no work has been performed investigating 
CaC >2 as a stable oxide of calcium at high pressures. Some 
previous work has explored the effect of low pressures 
(<10 GPa) on the bond lengths and lattice parameters 
of /4/mmm-Ca02 [5]. We therefore employ DFT calcu¬ 
lations to explore the behaviour of Ca 02 at pressures in 
the GPa range. DFT calculations provide an excellent 
avenue for investigating materials properties under pres¬ 
sure, both at pressures accessible to diamond anvil cells 
[8] and at terapascal pressures [9, 10]. To explore the 
stability of CaC> 2 , we search for new crystal structures 
of this compound at a variety of pressures in the range 
0-200 GPa. 


METHODS 

Density functional theory calculations are performed 
using the CASTEP plane-wave pseudopotential code 
[11]. Ultrasoft pseudopotentials [12] generated with the 
CASTEP code are used for both calcium and oxygen, 
with core states ls 2 2s 2 2p 6 and Is 2 , respectively. We 
use the Perdew-Burke-Ernzerhof (PBE) [13] form of the 
exchange-correlation functional with a plane-wave basis 
cutoff of 800 eV. A fc-point sampling density of 2-7T x 0.03 
A -1 is used for our CaO and Ca02 phases. For our 
oxygen phases, we use a denser /e-point sampling of 
27t x 0.02 A -1 . Bulk modulii are calculated by fitting 
static lattice pressure-volume data to the third-order 
Birch-Murnaghan equation of state. 
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The electronic density of states is calculated using the 
OPTADOS code [14-16]. Calculations of phonon frequen¬ 
cies are performed with the CASTEP code and a finite- 
displacement supercell method, using the quasiharmonic 
approximation [17, 18]. 

To search for new phases of CaC> 2 , we use the ab ini¬ 
tio random structure searching (AIRSS) technique [19]. 
AIRSS proceeds by generating random starting struc¬ 
tures containing a given number of formula units. Of 
these, structures which have lattice parameters giving 
reasonable bond lengths and cell volumes at a particu¬ 
lar pressure are then relaxed to an enthalpy minimum, 
and the lowest enthalpy structures are selected for refine¬ 
ment. AIRSS has proved to be a very powerful tool in 
predicting new structures, several of which have subse¬ 
quently been found experimentally. AIRSS searches have 
for example uncovered high pressure phases of silane [20] , 
and correctly predicted high pressure metallic phases in 
aluminium hydrides [21]. 

In this study, we perform structure searching at pres¬ 
sures of 0, 10, 20, 50, 100, 150 and 200 GPa. The bulk 
(about 60%) of our searches use cells with 2 or 4 formula 
units of CaC> 2 , and we have also performed searches with 
cells containing 1, 3, 5, 6 and 8 formula units. Not all 
combinations of pressures and formula unit numbers are 
searched. In total, we relax over 25,000 structures in 
our CaC >2 searches. We supplement our AIRSS searches 
by calculating the enthalpies of five known alkaline earth 
metal peroxide structures taken from the Inorganic Crys¬ 
tal Structure Database (ICSD) [22], and other authors 
[7, 23, 24]. The relevant alkaline earth metal is replaced 
with calcium where appropriate. 


CaO 

CaO undergoes a transition from the rocksalt ( Fm3m ) 
to the CsCl (Pto3to) structure (the B1-B2 transition) 
around 60-65 GPa. Diamond-anvil-ccll experiments in¬ 
dicate a transition pressure of 60±2 GPa at room temper¬ 
ature [25] , while DFT calculations give a transition pres¬ 
sure of 65-66 GPa [26, 27]. We find a pressure of 65 GPa 
in the present study, and we therefore use the Fm3m 
rocksalt structure below 65 GPa and the CsCl structure 
at higher pressures. The bulk modulus of Fm3m -CaO 
has been measured to be 104.9 GPa [28], while our static- 
lattice DFT calculations give a value of 108.5 GPa. To 
exclude the possibility that CaO might have a different, 
more stable structure (other than Fm3m or Pm3m) over 
the pressure range 0-200 GPa, we also perform AIRSS on 
CaO at pressures of 50, 100 and 200 GPa. We do not find 
any new low-enthalpy structures for CaO at these pres¬ 
sures. 


Solid oxygen 

Structure searching over the pressure range 0-200 GPa 
has already been performed for solid oxygen [29, 30], 
and we use the lowest enthalpy structures found therein. 
We also examine the enthalpies of the experimentally- 
determined a and S oxygen phases at low pressures [31]. 
Choosing the lowest-enthalpy oxygen structure at each 
pressure, we find that S- O 2 is stable between 0 and 1.2 
GPa, after which an insulating phase with space group 
Cmcm [31] is stable up to 41 GPa. A phase of symmetry 
C2/m [30] then becomes stable, remaining so up to 200 
GPa. Our spin-polarised calculations show no discern- 
able difference in enthalpy between a <5-02 phase with 
antiferromagnetic spin ordering, and the experimentally- 
determined ferromagnetic spin-ordering for 8 -O 2 [32]. 

These results are not in accord with low-temperature 
experiments on solid oxygen, which predict a -02 to be 
stable between 0 and about 5 GPa, <5-02 to be stable 
between about 5 and 10 GPa, followed by the ‘e-02 5 
phase between 10 and 96 GPa, with a further isostruc- 
tural phase transition around 100 GPa [33-35]. Our DFT 
calculations do not yield a- 02 ; optimising its structure at 
low pressures simply gives the 8-0 2 structure. The afore¬ 
mentioned C 2/to oxygen phase is however very similar in 
structure to e-02 (which is also of C2 /to symmetry), and 
is within 50 meV/f.u. of that phase in enthalpy over the 
pressure range 0-200 GPa. Any higher enthalpy structure 
for oxygen over the pressure range being explored here 
would only increase the calculated stability of Ca 02 , so 
we proceed with the lowest enthalpy DFT phases for oxy¬ 
gen. 


RESULTS 

Structure searching and static lattice results 

Our structure searching is carried out by minimising 
the enthalpy at a given pressure within the static-lattice 
approximation. We discuss these results first before pre¬ 
senting our calculations of the Gibbs free energy. 

Fig. 1 shows the enthalpy-pressure curves for nine 
phases of Ca02- The first five of these, labelled with 
their space group symmetries /4/tototo, Pa3, Pna2\, 
Cmmm and /4/tocto in Fig. 1, are known alkaline earth 
metal peroxide structures as mentioned in the ‘Methods’ 
section. The other four, with space group symmetries 
C2/c and P2\/c , are new Ca02 structures. These are the 
lowest-enthalpy phases that turned up during our AIRSS 
searches. The dotted line in Fig. 1 shows the enthalpy of 
CaO+^ 02 , calculated using the lowest-enthalpy phases 
of CaO and O 2 at each pressure. Any Ca02 phase below 
this dotted line is stable against decomposition in the 
manner of Eq. (1). 
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FIG. 1: Static lattice enthalpies, in eV per unit of Ca02, of calcium peroxide phases in the pressure range 0-5 GPa (left) and 
0-200 GPa (right). Enthalpies are given relative to the IA/mmm phase of Ca02- C2/c-I, C2/c-II, P2\/c-H and P2i/c-L are 
new structures of Ca02 found using AIRSS. The left-hand plot highlights the enthalpy differences between the Pna2\ and two 
C2/c phases at low pressures; C2/c-l is the lowest-enthalpy structure at 0 GPa from our searches. The arrow and box in the 
right-hand figure indicate the scope of the left-hand figure. Structures of Ca02 below the dotted line are thermodynamically 
stable against decomposition into CaO and O 2 . 


We find the following sequence of phase transitions 
for CaC >2 by considering the lowest enthalpy structure 
at each pressure: 


/ ^ 1r . , T 1.8 GPa D 0 2.0 GPa , TT 
(72/c-I->• Pna2 1 ->• C2/C-II 

27.6 GPa ,. , 37.9 GPa , r 

- > IA/mcm - P2i/c-L, 


with the transition pressures between each structure as 
indicated by the arrows. The P2i/c-L phase is then pre¬ 
dicted to be stable up to at least 200 GPa. 

At 0 GPa, our structure searching reveals a number of 
Ca02 phases which are very close (within 10 meV/f.u.) 
in enthalpy. Within the present approximations (DFT 
with PBE exchange-correlation), a phase of C2/c sym¬ 
metry (‘C2/c-P) is lowest in enthalpy at 0 GPa. This 
phase is 8.4 meV/f.u. lower in enthalpy than the pro¬ 
posed Pna2i-symmetry ground state structure of Zhao 
et al. [7], which we note also turned up in our AIRSS 
searches. Our searches also uncover a second C2/c- 
symmetry phase (‘C2/c-IP) which is slightly higher in 
enthalpy than Pna2i at 0 GPa. The enthalpies of these 
three phases at low pressures are shown in more detail in 
the left-hand panel of Fig. 1. 

We find that the calculated XRD patterns for the 
C 2/c-I, C2/c-II and Pna2\ structures share many similar 
features with available experimental data [23], although 
this XRD data is best fitted by the Pna2\ structure. The 
enthalpy differences between the C2/c-I, -II and Pna2\ 
phases at 0 GPa are dependent on the choice of exchange- 
correlation functional. In Table I, we show the enthalpies 
of the Pna2\ and C2/c-II phases relative to the C 2/c-I 


phase using a variety of different functionals. We find 
differences in calculated equilibrium volumes as well: for 
C2/c-l at 0 GPa, the LDA gives a volume of 35.8 A 3 /f.u., 
while the PBE functional gives 39.1 A 3 /f.u. The LDA 
typically overbinds in DFT calculations, with PBE in¬ 
stead underbinding, so these two volumes likely bracket 
the true volume of the C2 /c-I phase at 0 GPa. Given that 
dP/dV is about -2.3 GPa/A 3 for this phase at 0 GPa, 
this volume difference (due to functional choice) corre¬ 
sponds to a pressure uncertainty in the neighbourhood 
of ±3 GPa. In light of this uncertainty, C2/c-I,C'2/c-II 
and Pna2\ are all reasonable candidates for the structure 
of CaC >2 at 0 GPa. 


TABLE I: Enthalpies of the C2/c-ll and Pna2\ phases of 
Ca02 at 0 GPa, in meV per unit of Ca02, relative to the 
C2/c-l phase using different exchange-correlation functionals. 


Functional LDA PBE PBESOL PW91 

WC 

C2/c-I 

0.0 0.0 

0.0 

0.0 

0.0 

Pna2\ 

0.1 8.4 

2.9 

8.7 

3.5 

C2/c-II 

-16.3 12.2 

-4.1 

12.0 

-2.7 


The C2/c-I, -II and Pna2\ phases exhibit very similar 
structures. Looking down the &-axis of each phase, cal¬ 
cium atoms form a nearly-hexagonal motif, and we note 
that the monoclinic angle /? is close to 120° for C2/c-l 
and C2/c-II. For the (72/c-I and -II phases, the peroxide 
ion axes are almost coplanar, as can be seen in Fig. 2. In 
C 2 /c-II, the axes of peroxide ions in the same plane are 
parallel, while in C2/c-I, they alternate in orientation in 
the same plane. We use red and blue colouring in Fig. 2 
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FIG. 2: 2x2x2 slabs of the C2/c-I (left) and (72/c-II (right) 
structures, viewed almost down the 5-axis. Both structures 
are very similar, with Ca atoms forming an almost-hexagonal 
motif when viewed from this angle. Green atoms correspond 
to Ca atoms, while red and blue correspond to O atoms. All 
O atoms in peroxide ions with parallel 0-0 axes are given the 
same colour. 

to show peroxide ions with parallel axes. 

AIRSS searches also reveal a second phase for CaC >2 
with P2i/c symmetry (‘P2i/c-iL’). As can be seen in the 
right-hand panel of Fig. 1, our static-lattice calculations 
show that this is not a stable phase for CaC >2 over the 
pressure range 0-200 GPa. However around 38 GPa, the 
enthalpy of this phase lies within 10 meV/f.u. of the 
M/mcm and P2\/c-L phases, opening up the possibility 
this phase could become stable once we take temperature 
into account. We consider this further in a later section 
(‘Lattice dynamics’). The IA/mcm phase reported here 
is also predicted for MgC >2 above 53 GPa [24]. 

We provide lattice parameters and bulk modulii for the 
C2/c-l , C2/c-II, /4/racm, P2\/c-H and P2\/c-L phases 
at 0, 20, 30 and 50 GPa in Table II. 

Bonding and electronic structure 

The bandstructure of the P2i/c-L phase at 50 GPa 
is shown in Fig. 3. The insulating nature of this phase 
is evident, with a calculated thermal bandgap of 2.4 eV 
and optical bandgap of 2.5 eV. These will be underesti¬ 
mates of the true bandgap owing to the use of the PBE 
functional. The lowest set of 4 bands is comprised almost 
entirely of Ca 3s orbitals. Above this, we find a central 
peak flanked by two smaller sidepeaks in the total density 
of states. The central peak is built from 12 bands and 
is largely Ca 3 p orbitals, while the two sidepeaks con¬ 
tain 4 bands apiece and are dominated by O 2s orbitals. 
Thus, we have a splitting of the O 2s orbital energy lev¬ 
els, possibly arising from the covalent bonding present 
in the peroxide [O-O] ion. Above this, and just below 


the HOMO (highest occupied molecular orbital), are 20 
bands which arise largely from O 2 p orbitals. The first or¬ 
bitals above the Fermi level consist of (unoccupied) O 2 p 
orbitals, followed by a dense band of Ca 3d orbitals. An 
almost identical pattern of bonding and electronic den¬ 
sity of states are found in our other low-enthalpy phases. 



FIG. 3: Bandstructure and electronic density of states of the 
P2i/c-L phase at 50 GPa. The density of states is shown 
projected onto the s, p and d angular momentum channels. 
We calculate a thermal bandgap of 2.4 eV, and an optical 
bandgap of 2.5 eV. The Fermi level is shown as a black dashed 
line. 



FIG. 4: Phonon dispersion relations and density of states for 
P2i/c-L Ca02 at 50 GPa. 

Fig. 4 shows the corresponding phonon band struc¬ 
ture and density of states for P2i/c-L at 50 GPa. The 
lack of imaginary phonon frequencies indicates the sta¬ 
bility of this particular phase. The distinctly separate 
high frequency bands around 930 - 1000 cm -1 in Fig. 4 
correspond to phonon modes that stretch the 0-0 co¬ 
valent bond in the peroxide ions. Two distinct perox¬ 
ide bond lengths are found in the ions of this structure: 
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TABLE II: Structures of the C2/c-I, (72/cTI, 14/mon, P2\/c-H and P2\/c-L phases of CaCH- 


Pressure 

(GPa) 

Space group 

Lattice parameters 
(A, deg.) 

Atom 

Atomic coordinates 

x y z 

Wyckoff 

site 

Bo (GPa) 

0 

C2/c (#15)“ (I) 

a=7.041 

6=3.685 

c=6.820 

Ca 

0.0000 

0.6399 

0.2500 

4e 

87 



a=90.0 

/3=117.8 

7=90.0 

O 

0.2548 

0.1404 

0.4119 

8f 


20 

C2/c (#15)“ (II) 

a=6.829 

6=3.403 

c=6.407 

Ca 

0.0000 

0.3413 

0.2500 

4e 

89 



a=90.0 

/3=118.8 

7=90.0 

O 

0.1661 

0.1553 

0.0252 

8f 


30 

M/mcm (#140) 

a=4.521 

6=4.521 

c=5.745 

Ca 

0.0000 

0.0000 

0.2500 

4a 

114 



a=90.0 

/3=90.0 

7=90.0 

O 

0.1143 

0.6143 

0.0000 

8h 


30 

P2 1 /c-H (#14) 6 

a=6.590 

6=4.842 

c=3.795 

Ca 

0.0611 

0.7628 

0.2781 

4e 

93 



a=90.0 

/3=105.2 

7=90.0 

O 

0.1301 

0.2714 

0.3087 

4e 







O 

0.2446 

0.4378 

0.1012 

4e 


50 

P2 1 /c-L (#14) 6 

a=4.223 

6=4.279 

c=5.949 

Ca 

0.0613 

0.5133 

0.2539 

4e 

no 



a=90.0 

/3=98.7 

7=90.0 

O 

0.0820 

0.0235 

0.4019 

4e 







O 

0.1224 

0.1182 

0.0089 

4e 



a C'12/cl - International Tables, Volume A: unique axis b. cell choice 1. 
b P12i/al - International Tables, Volume A: unique axis 6, cell choice 3. 


1.44 A and 1.46 A, which splits these higher frequency 
bands. These peroxide 0-0 bond lengths are somewhat 
longer those found in molecular oxygen, which has an O- 
O bond length of 1.207 A [23] at ambient conditions, and 
the bond length in C2/m oxygen at 50 GPa, which we 
calculate as 1.20 A. These longer 0-0 bond lengths are 
however typical of crystalline ionic peroxides [5]. 

Stability of Ca02 

As can be seen in Fig. 1, our low-enthalpy phases for 
Ca 02 show remarkable stability against decomposition as 
pressure increases. We predict Ca02 to be most stable 
at the B1-B2 CaO phase transition pressure of 65 GPa. 
There, the decomposition enthalpy (Eq. (1)) of Ca02 is 
+0.64 eV/unit of Ca02 (+62 kJ/mol). 

One implication of the stability of Ca02 is that it may 
be preferentially formed over CaO in an oxygen-rich envi¬ 
ronment under pressure, through the reverse of Eq. (1). 
For example, the pressure in the Earth’s lower mantle, 
at a depth of around 1550 km, is about 65 GPa [25], 
close to our predicted peak stability pressure for Ca 02 - 
The formation of Mg02 in this way has also been dis¬ 
cussed [24], although much higher pressures (>116 GPa) 
are needed before Mg 02 is stable against decomposition, 
whereas Ca02 is stable from around 0 GPa. The mantle 
temperature in the Earth at 65 GPa is in the neighbour¬ 
hood of 2500 K, and our phonon calculations show that 
under these conditions, AG = +0.54 eV/f.u. for the re¬ 
action of Eq. (1). Hence, Ca02 is a thermodynamically 
stable oxide at temperatures and pressures encountered 
in planetary interiors. 

Reactions of the form X+O 2 —> Y for species X and 
Y, such as the reverse of Eq. (1), are known as redox 
buffers. Such reactions are key in determining planetary 
mantle compositions. In the Earth’s mantle, there are 


a number of such buffers, usually involving the further 
oxidation of iron and nickel compounds, such as FeaCH 
+ JO 2 —> §Fe 203 and Ni + |02 —> NiO. The natural 
formation of Ca02 by further oxidation of CaO in Earth’s 
mantle, while energetically favourable at high pressures 
and temperatures, would also need to compete against 
the further oxidation of other such compounds. The av¬ 
erage CaO content in the Earth’s mantle is about 3% 
[28], though exoplanet mantles offer a rich variety of al¬ 
ternative compositions. 

We highlight the fact that the pressures at which these 
different Ca 02 phases are predicted to become stable are 
amenable to experimental study in diamond anvil cells. 
Very few pressure-induced phase transitions for A[B 2 ] 
compounds are experimentally known [36], although at 
least one example already occurs among the alkaline 
earth metal peroxides, in Ba02 [36]. Equivalently, it 
would be interesting to test the reactivity of CaO and O 2 
under conditions of excess oxygen, as our results suggest 
that CaO and O 2 are reactive at high pressures. The for¬ 
mation of Ca 02 may require laser heating in a diamond 
anvil cell to overcome the likely high potential barriers 
between phases. High pressure phases of Ca02 could be 
recoverable at lower pressures, although possibly not at 
ambient or zero pressure. 

Lattice dynamics 

In addition to our static-lattice calculations, we calcu¬ 
late the phonon free energies of our lowest-enthalpy Ca 02 
phases, namely those with Pna2i, C2/c-I/II, /4/mcm, 
P2\/c-H and P2\/c-L symmetries. We encounter no 
imaginary phonon frequencies over the pressure ranges 
relevant to these phases, indicating that they are dynam¬ 
ically stable. The relevant thermodynamic potential is 
now the Gibbs free energy, which includes the phonon 
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FIG. 5: T-P phase diagram of Ca0 2 as calculated using the 
quasiharmonic approximation and our lowest-enthalpy struc¬ 
tures. 


pressure. Selecting the lowest Gibbs free energy struc¬ 
ture from these six structures gives rise to the T-P phase 
diagram given in Fig. 5. We note that the upper-left 
(low-P, high-T) part of the phase diagram is representa¬ 
tive only, because at ambient pressures Ca0 2 decomposes 
at temperatures around 620 K [2, 3]. 

We find that the small enthalpy difference between the 
IA/mcm and P2i/c-H phases seen in our static-lattice 
calculations (Fig. 1) closes with increasing temperature, 
and we see the emergence of P2\/c-H as a stable phase 
for CaC >2 at 37.7 GPa and for T > 281 K. The free en¬ 
ergy difference between the IA/mcm and P2\/c-H phases 
does remain quite small however, around 20 meV per 
Ca0 2 unit at the most over the temperature range 0- 
1000 K. 

At room temperature (300 K), we therefore predict 
a different sequence of phase transitions than those for 
our static-lattice calculations. We find that, with PBE 
exchange-correlation: 


, T 2.2 GPa D 0 2.6 GPa , TT 
C2/c-I - Pna2\ ->• (72/c-II 

r , , 37.4 GPa „„ , „ 37.9 GPa 

M/mcm - > P2\/c-H ->• 


28.3 GPa 
- > 

P2i/c-L, 


with the arrows labelled by the predicted transition pres¬ 
sures. 

The phase diagram of Fig. 5 does not extend all the 
way to 200 GPa. However, we expect P2i/c-L to con¬ 
tinue to be the most stable phase at high pressures. 
This is because our structure searching results (which use 
static-lattice enthalpies) show that the next most stable 
structure for Ca02 over the pressure range 100-200 GPa 


is at least 45 meV per unit of CaC >2 higher in enthalpy. 
This was not the case at low pressures, where our searches 
reveal quite a few structures (such as P2\/c-H) that are 
close to becoming stable and may therefore do so at high 
temperatures. 


CONCLUSIONS 

Structural changes in Ca02 under pressure have been 
explored over the pressure range 0-200 GPa at tempera¬ 
tures up to 1000 K. CaC >2 remains insulating up to pres¬ 
sures of at least 200 GPa. Structure searching and DFT 
calculations reveal six stable phases for Ca 02 over these 
pressure and temperature ranges, of which five are re¬ 
ported for the first time in this study. Calculations of the 
phonon frequencies of these new structures confirms their 
dynamical stability. The lowest-enthalpy phase of Ca02 
at 0 GPa is dependent on the choice of DFT exchange- 
correlation functional. At pressures above 40 GPa, a 
phase of P2 1 /c symmetry ( l P2i/c-L') emerges for Ca0 2 
which is predicted to be stable up to 200 GPa, and at 
mantle pressures and temperatures. Ca 02 is a very sta¬ 
ble oxide of calcium at high pressures, and may be a con¬ 
stituent of exoplanet mantles. The pressures at which 
these Ca 02 phases become stable are readily attainable 
in diamond anvil cells. 
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